Ultrasmall double junction in terms of orthogonal polynomials 
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The "orthodox theory" of a single electron double junction is dealt with. It is shown that 
the stationary solution of the underlying master equation allows the construction of any time- 
dependent solution in terms of orthogonal polynomials. The approach pays off if the stationary 
solution becomes simple. Two special cases are considered. We use the time-dependent solution to 
calculate the current noise in these cases. 

73.40Gk, 73.40Rw 



I. INTRODUCTION 

Due to its simplicity the normal metal single electron 
double junction (1E2J, shown schematically in Fig. 3) is 
one of the best-studied systems of single electronics, both 
experimentally and theoretically. The.,common theoreti- 
cal approach uses "orthodox theory" 00 to find a master 
equation for the charge probability of the island between 
the two junctions. The stationary, solution of this master 
equation is known analytically.oEl Since the derivation of 
the equation itself restricts its solution to slow processes 
on the time scale of i?sCs = (Ri +R%){C\ +C2) the sys- 
tem appears more or less transparent within "orthodox 
theory." 

However, if one applies "orthodox theory" to the 1E2J 
one realizes that the structure of the stationary solu- 
tion restricts its practical use to both low voltages and 
temperatures where the energy scale is the Coulomb en- 
ergy Ec = e 2 /(2Cs). Therefore the investigation of n- 
level systems, taking into account-the most important 
charge states only, was introduced!] Using this method, 
main features of single electron tunneling (SET), such 
as Coulomb blockade or the Coulomb staircase, can be 
understood theoretically. 

Here we present an alternative method of solving the 
master equation that is mainly applicable at either high 
voltage and negligible temperature or high temperature 
and small voltage. In these cases the stationary solution 
of the master equation simplifies and provides access of 
the full time-dependent solution. In this regard the pre- 
sented approach is complementary to that of a n-level 
system. We use this solution for the calculation of the 
current and the current noise in a 1E2J for frequencies 
/«(i? E C s )- 1 . 

The physical relevance of our results lies in clarifying 
the asymptotic behavior of the 1E2J rather than in un- 
veiling new SET effects. In fact, the latter can not be 



expected in the considered domain. Furthermore, the 
treatment may be used as an easy-to-handle approxi- 
mation in specific cases. 

In Sec. H we present the general approach. In Sees. Ill 
and |y] this general approach is applied to the case of a 
high voltage at low temperature (ohmic limit) and low 
voltage but high temperature (thermal limit), respec- 
tively. Both stationary mean current through a 1E2J 
and current noise power spectrum are calculated. The re- 
sults are discussed and limiting cases are compared with 
known results. 



II. METHOD 

Let r M (x) and Z M (x) denote the tunneling rates in right 
and left direction (see Fig. |l|) across the /i-th junction if x 
excess charges reside on the island between the junctions. 
Then the matrix 



W(y, x) = (-l) 1+x +y{ [ n (x) + l 2 (x)] [S yx+1 + 5 
+ [r2(x) + h{x)] [8 yx + S yx+ i] } 
enables to write the "orthodox" master equation 
\&)=W\a), 



(1) 



with the probability cr{x) of x excess charges on the is- 
land as the x-th component of the vector |cr). Due to 
the underlying tunneling of charges x is integer. The left 
and right eigenfunctions are given by 

W|<£ n ) = -\ n \cj> n ), 
(V'„|W = -A„(^|, 

with the eigenvalues A„ > 0. It can be showni that one 
of the eigenvalues equals zero (let it be Ao). The corre- 
sponding cigenfunction |<^>o) describes the stationary so- 
lution of the master equation ([!]). It obeys the detailed 
balance condition 
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W(ar, x + l)Mx + 1) = W(a! + l t x)M*) (2) 
Left and right eigenfunctions are related throughi 

4> n (x) = 1p n (x)(f> (x). (3) 

From Eq. (||) we see that ij)o(x) = 1. The orthogonality 
relation of the left eigenfunctions reads 



^i/j m (x)^„(x)(f) (x) = S n 



(4) 



Thus, the eigenfunctions ipn(x) form a set of orthogonal 
polynomials with the weight function <po(x) ("SET poly- 
nomials"). As an example, these polynomials are shown 
in Table | for a special case. In general, ip n (x) is a poly- 
nomial of order n in x. 

If the stationary solution of the master equation (Q) is 
known, an arbitrary time dependent one follows readily 
from Eqs. (§) and (§: 



\a(t))=J2b n e 



X n (t-t ) 



where the b n are adjusted to meet the initial probabil- 
ity distribution \a(to)), b n = (ipnW(to)) ■ The eigenvalues 
A n , which determine the time-dependence, result from 
Eq. (Eh: 



A„ = W(x-l,a:) 
+W(x + l,x 



1 



1p n (x - 1) 
1pn(x) 

_ 4> n (x + 1) 

1pn{x) 



(5) 



Vx. 



In summary, the stationary solution of the master equa- 
tion complemented with an initial condition determines 
the full time-dependent solution. 



III. OHMIC LIMIT 

Now we turn to the case of zero temperature. Then 
the detailed balance property (|j) causes the ratio 4>o (x + 
l)/4>o(x) to be a rational function of the bias voltage v = 
V Cs/e, and of the dimensionless charges q = Q/e and x. 
Whereas x describes the integer part of the charge, which 
is altered during tunneling processes, the parameter q 
expresses the influence of background charges and/or an 
applied gate voltage. In the following we will make use of 
the additional dimensionless parameters k m = Ca—^/Cs 
and = i? At /i?s, [i — 1, 2. If v is integer and obeys 



+ q — 1/2 = to — 0^ 



(6) 



for an arbitrary other integer to, the ratio <f>o(x+l) / <f>o(x) 
can be written as 

M x + 1 ) 62 N - k 



<j) (x) 



Qt k + 1 ' 



with N — v — 1 and k — x — x min — x + m — 1. 
^min denotes the lowest of all accessible charge states, 
a^min = —[k-2V + q + 1/2], and [. . .] the integer part of . . .. 
Therefore 4>q(x) is a binomial distribution in these cases: 



(j) (x) 



N 



The given condition (^) defines a number of points in 
the V-Q plane, as shown in Fig. || ("binomial points"). 
However, if 1 k <C N, the condition (||) is fulfilled for 
all v and q. This is the case for high bias voltage V. 

For the binomial distribution the polynomials ip n {x) 
are known. They caa be expressed in terms of 
Krawtchouk polynomials^] as follows: 



ip n (x) = K n (x ~ x min ; Q 2 , N) 



The eigenvalues are determined from Eq. (^|), which is 
the difference equation (1.10.4) of Ref. 0, yielding 



n 1 



Ry>C-£, Q1Q2 -RyCs 



(7) 



with #|| = R^/Rt,- 

The calculation afj-,the stationary current follows the 
standard procedureH'S and yields 



(I) 



e(v - 1) 
Ry,Cy, 



(8) 



One verifies that this formula describes the asymptotic 
behavior for large bias correctly. For lower bias the "bi- 
nomial points" of Fig. || correspond to voltages where the 
(I)-V characteristic of the 1E2J touches the asymptotic 
line. For the symmetric 1E2J with q = this happens 
for v = 2 m + 1 with integer m. 

Next we present our result of the current noise of a 
1E2J. The method of calculation has been described in 
Refs. [l(] and [ll]. It takes advantage of the representation 
of the conditional probability p(m, t; to', t') and the ma- 
trix elements of W in terms of the eigenfunctions ip n (m). 
After several steps (see Appendix) we obtain 



S(0) = 2e(I)(l-2 Ql g 2 ) 



(9) 



for zero frequency. Since a realistic estimation shows that 
ORsCs)- 1 > 1GHz typically, the result (|) is valid in a 
wide range of frequencies, down to the regime where 1//- 
noise becomes predominant. 

Eq. (^]) yields the ordinary shot noise formula in the 
limit of an extremely asymmetric 1E2J (pi — > 1 A Qi — > 
or vice versa), where the current across both junctions is 
uncorrelated. The shot noise is halve d for , a symmetric 
1E2J in agreement with earlier results.113113 

Additional noise suppression in double junctions in 
dependenee|-ep the bias voltage was both ,-tiieoretically 
predicted3li2l and experimentally observedJi-2 Our result 
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(0) does not recover this suppression since the "binomial 
points" are located at bias values where this suppression 
does not occur. 

In Fig. | and Fig. ^ the results for symmetric and asym- 
metric 1E2J are shown. For the symmetric system the 
noise is suppressed. An additional suppression occurs 
between the "binomial points" in comparison to 2e(I). 
For higher bias the use of (^|) as approximation improves. 
This improvement is best for the symmetric 1E2J. 



IV. THERMAL LIMIT 

Let us now consider the case of dominating tempe- 
rature k^T = > Eq. The equilibrium distribu- 
tion is then a Gaussian with mean (/ti — Q\)v and width 

(2/3£ C r i/2 



(j>o{x) 



I3E C 



■ exp i 



[x-( Kl - gi )v] 2 pE c ), (10) 



which follows from detailed balance (g), or equivalently, 
from a thermodynamic consideration. In order to de- 
termine the eigenvalues, we consider Eq. (||) in the limit 
x — > ±00, where it yields again the temperature indepen- 
dent values of Eq. ([?]). While our consideration has been 
exact within orthodox theory so far, the treatment re- 
quires several approximations from now on. Since many 
states are significantly occupied, the above charge state 
sums can be replaced by integrals J2 X ~~ * I Then, 
the SET polynomials are given in terms of Hermite poly- 
nomials 

ip n (x) = i == '-. 11) 

V2 n n\ 

For the calculation of the stationary current let us in- 
troduce a current vector (t| by 

(J) = (,1a) 
2 



where a Taylor expansion in (/3-Ec) is used to approxi- 
mate the high temperature behavior of (t|. Even if this 
expansion yields positive powers in x, which are inte- 
grated later on, its use is justified by <f>o(x) of (|Tc|), which 
decreases very rapidly with increasing |x|. For the sta- 
tionary current we find 



(I) = (i\<f>o) = a Q 



*o V_ 



1 



(3E C 



This high temperature result is independent of the dou- 
ble junction's symmetry|— It agrees with the known result 
from SET thermometry]^ 



The calculation of the noise follows again Refs. 0Jll 
In terms of rti 



M x±1 ) = M x )e~ 0Ec{1±2lx ' (Kl ' ei)v]} - 

Furthermore, high temperature expansions are used 
again to simplify the calculation which is lengthy, but 
straightforward. Restricting to the leading term and its 
first order correction results in contributions from the 
eigenvalues A .i only. After carefully collecting all rele- 
vant contributions the zero frequency noise power 



3(0) = 



4fc B T 



(3E C 



2- 



Ql 



02 



(12) 



is found. This result coincides with the high tempera- 
ture limit of Ref. |l^, but gives a correction to the known 
expression as well. In case of symmetric 1E2J (or more 
generally: for Ki — Qi) the term (...) in ( |l2] ) simplifies 
to unity and it becomes clear that the Nyquist formula 
yields an upper bound of the noise of a 1E2J. This fact 
expresses the partial coherence between tunneling events 
on both junctions. At high temperature, however, this 
coherence get lost— Hiis behavior does not correspond to 
the coth-formulaliirE 



5(0) = 2e(J) coth 



/ eV 
\2k B T 



using the current / through the system and the bias volt- 
age V, which predicts elevated noise above AksT/R^- We 
want to point out that the coth-formula was derived for 
the single junction (1E1J) and does not allow for coher- 
ence effects of both junctions. 



V. CONCLUSIONS 



We have presented a method allowing access to the 
time-dependent solution of the master equation that is 
used in the semiclassical "orthodox theory" to describe 
the behavior of an ultrasmall single electron double junc- 
tion. This method is based on general properties of mas- 
ter equations. We can show that the stationary solu- 
tion of the master equation is sufficient to determine any 
time-dependent solution by quadratures on the range of 
allowed charge states. For special cases of the stationary 
solution this procedure simplifies to an approximatively 
analytically solvable problem. Two of these cases are 
presented in this paper in detail. After calculating the 
respective mean currents the time-dependent solution is 
used to find expressions for the current noise power spec- 
tra. The obtained results are discussed. It turns out 
that the presented method together with the usual n- 
level approach enable access to time-dependent solution 
in a wide parameter range of the bias voltage V and the 
temperature T. 
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The eigenvalues (0) follow directly from the difference 
equation (1.10.4) in Ref. ^|. In the evaluation of the noise 
we follow strictly the description of Ref. [n], namely 



c-Sh _ 



^ 2 E 



+ ; M (n)]0 o (n), 



(Al) 



S c ^(t-t')=e 2 J2 [rM)-Urn)]p{rn,t-,m',t') (A2) 

m m' 

[n(m' - l)4>o{m' - 1) - h{m! + l)0 o (m' + 1) + r 2 (m' + l)4>o{m' + 1) - h{m! - 1)0 o (to' - 1)] , 



S^(t) = -(I) 2 + 5(t)5^S* h + S^,(t), (A3) 
S(t) = J2 K ^' s w'( t )> ( A4 ) 



S(lo) = 4 / dtcos(ujt)S(t) 
Jo 



(A5) 



is finally the noise power spectrum of the current that is 
considered in the zero frequency limit u) — ► 0. 

The conditional probability p{m, t] to' , £') can be ex- 
pressed in terms of the eigenfunctions of W 



JY 



p(m, t; to', t') = 2J i(Jn(m')(t) n (m)e~ 



-A„(t-t') 



From the orthogonality relational 

^ ip n (rn')<p n (m) = 5 mm > 

n=0 

we see that the correct initial condition of p{m, t; to', i') 
is obeyed at t = i'. Now the different noise contribu- 
tions are calculated. Due to the vanishing temperature 
one finds S^ h = e(I) independent of the junction num- 
ber \i and the system's symmetry (/cj 2i £1,2)- The non- 
diagonal terms of S^, simplify using the detailed balance 
property (||) . With the rates ry, (to) (the current direction 
has to be considered only) Sf . , is transformed into 



N 



^<(i-i')=e 2 ]Te- 



\„(t-t') 



r M (m)(f>n(m)r3- M >(m')(l>n(rn'). 



The term corresponding to n — in the last equation re- 
sults in (/) 2 , whi ch c ancels against the same term in the 
noise expression (A3). The rates r^(m) are linear func- 
tions of their arguments. Hence, they can be expanded 
into the left eigenfunctions tp n (m) with contributions of 



"00,1 (to) only. Note, that the 0-function of the rates can 
be omitted in the considered "binomial points" . Using 
the representation 77,(771) = c fJ _QipQ(m) + c /j iV'i(to) the 
expression above is further simplified to 

S"> (t - t') = e 2 c M lC3 _ M ; 1 exp[-Ai(i - t% 



In detail, we find 



.o = (I)/e, 



Cul 



(-l)H-l 



N 



Q1Q2 



using the self-duality property of the Krawtchouk poly- 
nomials. Finally, assemblage of the noise follows the line 
of Ref. [ll] straightforwardly. For zero frequency we ob- 
tain 



S(0) = 2e<J)(l-2eift,) = 



2e 2 (v- 1) 



(1-2^02), 



-Reg's 

where the stationary current is expressed via (| 



APPENDIX B: CALCULATION FOR THE 
THERMAL LIMIT 

For the calculation of the eigenvalues A„ the approxi- 
mations 



W(x - 1, as) = r 2 (x) + h(x) 



1 x 



Ry,Cz Q1Q2 
_ iJj„(x - 1) ^ n 

W(x + l,x) = r^x) +h{x) w 

are found for x ^> 1. Introducing into (jfy results in 
the eigenvalues ([?]). A similar consideration for x <C — 1 
yields the same result. 

The calculation of the mean current starts from the 
Gaussian approximation of the stationary solution 4>q (x) 
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of (10) and the corresponding left eigenfunctions (11). 
The relevant expansion of the current operator i(x, v) is 
simplified by the use of the shifted charge y — [x — [k\ — 
qi)v\: 



l{x, v) 



v_ 



e(«i ~ Q\)y 



(3E C 



The determination of <Zj results in contributions to ao,i, 
in detail 



«i 



= V (i P Ec \ 
e(/ti - gi) 



1 



The calculation of the stationary current makes use of ao 
only. 

The calculation of the noise follows a similar, though 
more complicated, line of the curr ent ca lcul ation above 
in the evaluation of the expression (Al)-(A5). By use of 
© we find 

ri,2(z T l)0o (a; T 1) = h,2(x)<t>o(x) exp(2/3E c gi^v) 
£1,2(2 ± l)4>o{x ± 1) = n t 2(x)(f>o(x) e-xp(~2f3E C Qi,2v). 

Furthermore, the abbreviations 

V(^) = e[r M (a;) - l^x)} 

3n{x) = e[l^(x) exp(2(3E c g ll v) - r^(x) exp(-2 ( 5£' c ^u)] 
k^(x) = e[r^(x) + l^x)], 

allow the transformation of the correlation function into 
S c ^(t-t') = X>P[-A„(i-*')]<vl0n) OVK) (B2) 
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and the Shottky value is expressed as S^ h — e(k fJi \(f> ) . 
The representation of the relevant high temperature ex- 
pansion of i^, j^, and k^ in terms of the left eigenfunc- 
tions ipn reads 



(vl 
O'mI 



2E C I, i 2i? c 
2E C 



1 



2 



ei? s Qf+VWEc] 
2E C 1 



1 



E c \, , . E C V2 



3eRT,g, 



0Ec 
3 

0E C 
3 



The calc ulati o n of the zero frequency noise power follows 
directly (A3)-(A5). Keeping the leading two orders in 
0Ec only results in the expression (E2I) given in the text. 



V/2 

o- 



R1 Ci Q ^2 P 



l 2 ^2 . V /2 



FIG. 1. Circuit diagram of the ultrasmall double-tunnel 
junction, consisting of two tunnel junctions connected in series 
(with resistances Ri,2 and capacitances Ci,2, respectively) 
1 hom@phys.unit.no. and an excess charge Q on the island. 
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FIG. 2. V-Q plane for C2/C1 = 1.5. The tiles correspond 
to regions of fixed N and the black corners indicate the points 
where the stationary charge distribution is binomial. 
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FIG. 3. Noise in a symmetric 1E2J (i? E = 200kfi, 
Ce = 0.2fF). The full line represents Eq. @ and the dashed 
line displays the 2e(I). The arrows indicate the "binomial 
points" for this system. 




V (mV) 

FIG. 4. Noise in an asymmetric («i = Qi = 0.9) 1E2J 
(Rv = 200kfi, C E = 0.2fF). The full line represents Eq. (|) 
and the dashed line displays the 2e(7). The arrows indicate 
the "binomial points" for this system. 



c„,o 



Cnl 



Cn2 



n = 
n = 1 

n = 2 



3c 



v + 1 



2r 




3u - 1 

V2v 2 - i 



TABLE I. Coefficients c„ m of the SET polynomials for a 
symmetric 1E2J with Q — T — at dimensionless voltages v 
in units of e/Cs so that v G (1;3]. The polynomials are built 
from these coefficients via ip n {x) = ^2 m c nm x m . 
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